Bound entanglement in the XY model 
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We study the multi-spin entanglement for the ID anisotropic XY model concentrating on the 
simplest case of three-spin entanglement. As compared to the pairwise entanglement, three-party 
quantum correlations have a longer range and they are more robust on increasing the temperature. 
We find regions of the phase diagram of the system where bound entanglement occurs, both at zero 
and finite temperature. Bound entanglement in the ground state can be obtained by tuning the 
magnetic field. Thermal bound entanglement emerges naturally due to the effect of temperature on 
the free ground state entanglement. 



Entanglement is a resource in quantum information 
science More recently it has become clear that un- 
derstanding the nature of quantum correlations may also 
help in a deeper description of complex many-body sys- 
tems (see Ref. [1] for a recent review on this topic). The 
work developed in these years in between the two areas of 
quantum statistical mechanics and quantum information 
has lead to several interesting results in both disciplines. 
Among the different aspects investigated so far we men- 
tion the study of entanglement close to a quantum phase 
transition Most of the work so far was developed 

to bipartite entanglement with some notable exceptions. 
On the other hand the entanglement monogamy prop- 
erty 0] constrains the entanglement sharing, and put 
an upper bound to the pairwise entanglement. 

Several indications demonstrate that the multipartite 
entanglement is indeed particularly important for the col- 
lective behavior of the system. It was shown, for example, 
that multipartite entanglement is enhanced with respect 
to the two-particle entanglement near a quantum crit- 
ical point 0, d, Q. Multipartite entanglement close to 
quantum phase transitions was quantified by the global- 
entanglement measure of Meyer and Wallach in [lO| or 



the geometric measure of entanglement The gener- 
alized entanglement measure introduced in [IJ] was used 
to study a number of critical and disordered 14 1 spin 
models. 

A comprehensive classification of the type of multipar- 
tite entanglement in spin systems has been recently given 
by [is,, 16]. Different bounds obtained to the ground state 
energy were obtained for different types of n-particle 
quantum correlated states. A violation of these bounds 
implies the presence of multipartite entanglement in the 
system. As discussed in [17| the analysis of the aver- 
age measures of multipartite entanglement might not be 
sufficient and the analysis of the distribution of block en- 
tanglement for different partitions may give additional 
information. 

In the present paper we focus on the Multiparticle En- 
tanglement (ME) with the aim to shed light on how en- 
tanglement is shared in a many body system. Specifically 



we analyse the type of entanglement of a subsystem made 
of few particles, tracing out the rest of the system and we 
study the ME shared between them (this approach was 
recently carried out also for the free electron gas jl8l.il9|). 
By this approach, although we cannot discuss the global 
ME properties of the whole system, we can gain insight in 
the details of such few-particle ME. In particular we con- 
sider the ID XY model in transverse field and we study 
the simplest multiparticle case of a subsystem made of 
three arbitrary spins of the chain. We analyze bipartite 
entanglement between a spin and the other two with re- 
spect to all possible bipartitions. We demonstrate that 
ME extends over a longer range than two-particle en- 
tanglement. An important feature, emerging from our 
analysis of ME, is the existence of Bound Entanglement 
(BE) 0,1211 shared among the spins of the chain. This 
peculiar form of entanglement is characterized by the im- 
possibility of distilling it into a pure form. It is a weak 
form of entanglement having features of both quantum 
and classical correlations. Bound entanglement has been 
a subject of intense research in the last years since it was 
shown to be a useful resource in the context of quantum 
information [H, [H 0, S [111 . 



Here we show that it emerges from the equilibrium 
properties of spin chains. Indeed for the ground state of 
quantum XY model spins far apart enough may be in 
a bound entangled state, even in the thermodynamical 
limit. At T ^ we find that ME is more robust than 
two-particle entanglement and increasing the tempera- 
ture it always turns into a thermal multiparticle bound 
entanglement. 



The paper is organized as follow: in the following sec- 
tion we review the scenario of three-qubit entanglement; 
we then briefly describe the model studied and some 
known results concerning its entanglement properties. In 
the last section we show our results at T = and at finite 
temperature. 
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FREE AND BOUND THREE-QUBIT 
ENTANGLEMENT 



For two-qubit systems, Concurrence [27| is a measure 
of entanglement that can be easily calculated both for 
pure and mixed states and it can be related to the 'en- 
tanglement of formation' [2^. Besides, any two-qubit 
entangled state can be converted into Bell states with 



the process of distillation [28|, 

For the entanglement of a three-qubit system, a much 
more complex scenario emerges. The total entanglement 
shared between the three qubits cannot be described just 
by studying the qubit/qubit entanglement between all 
the pairs, as measured f.i. by the Concurrence of the 
reduced 2-qubit density matrices. In fact in general the 
three qubits share a ME whose complete information is 
unavoidably lost by tracing out a qubit. The well known 
paradigmatic example are the GHZ states for which each 
couple of spins shares no entanglement, nonetheless each 
spin is maximally entangled with the other two [7|. A 
further qualitative difference that emerges for the three 
qubits entanglement is the existence of bound entangled 
states. The classification 'free' and 'bound' entangle- 
ment can be drawn with respect to distillation properties 
of the entanglement [2^ . Entanglement is bound (not- 
distillable) if no maximally entangled states between the 
parties of the system can be obtained with local oper- 
ations and classical communication (LOCC), not even 
with an asymptotically infinite supply of copies of the 
state 2l|. Despite thepragmatic definition of BE, the 
nature of its correlations is peculiar. In fact BE is a very 
weak form of entanglement that has both quantum and 
classical features. For instance some examples of multi- 
partite bound entangled state violating Bell-type inequal- 
ities were found [30|, but recently it was shown that if 
we allow collective manipulations, and postselection, no 
bound entangled state violates Bell-type inequalities 

While for two-qubit system no BE states exist 



29|, 



31|. 



the case of a three-qubit system the relative Hilbert space 
is large enough to allow such structure to appear. An in- 
teresting case in which such BE appears is related to the 
'incomplete separability' of the state (see 2]J). This con- 
dition happens for example when a state of a tripartite 
system A-B-C is separable with respect to the partition 
A — BC and B — AC and non-separable with respect to 
C — AB. The 'incomplete separability' is a sufficient con- 
dition for a state to have BE since the three qubits are 
entangled and no maximally entangled state can be cre- 
ated between any of the parties by LOCC. For exam- 
ple, no entanglement can be distilled between C and A 
because no entanglement can be created with respect to 
the partition A — BC by LOCC. In the following sections, 
the feature of incomplete separability will be exploited to 
detect such a kind of BE in the system. 

In order to analyze the free and bound entanglement 
between the qubits we focus on the entanglement between 



a qubit and the other two with respect to all possible 
bipartitions. To measure such bipartite entanglement we 
used the Negativity [s^ : 



(1) 



where are the negative eigenvalues of the partial trans- 
pose of the density matrix with respect to the bipartition 
v. Negativity gives the degree of violation of Peres con- 
dition and was proved to be a measure of entanglement 
[33 , [33! , though it can not detect Positive Partial Trans- 
pose (PPT) entangled states. 

We remark that the genuinely multipartite entangle- 
ment (e.g. GHZ type of entanglement) is distinct from 
generic multiparticle entanglement (e.g. a mixture of 
Bell states shared between the parties) jsj] • In (ssf a 
Negativity-based multipartite entanglement measure was 
proposed, but it is not easy to handle for generic mixed 
states since it requires a convex roof minimization. For 
a multiparticle system, Negativity is able to measure the 
'global' ME with respect to a certain bi-partition. For 
instance. Negativity can not distinguish between qubit- 
qubit entanglement and ME. For example, if qubit A is 
entangled separately with qubit B (as measured by a non 
zero concurrence Cab) then A will surely be entangled 
with B and C as a whole and thus Ma\bc 0- O'^ the 
other hand if no qubit/qubit entanglement is present in 
the system, then a non-zero Negativity detects multipar- 
ticle entanglement. Besides, carrying such analysis with 
respect to all possible bipartitions can give information 
on how such ME is shared between the qubits. It is pre- 
cisely this method that we will employ to study the ME 
entanglement in the system. 



THE XY MODEL 

Entanglement has been studied in a variety of spin 
models 2]. A lot of attention has been devoted to the 
ID quantum XY model. The Hamiltonian is defined as 

N 

H^JJ2{1 + l)SfSf+, -f (1 - 7)^r^^i - BSt. (2) 



The spins S*" = ^ <^ — ^^Ui ^ (o"" ^-re Pauli matrices) 
experience an exchange interaction with coupling J > 
and uniform magnetic field of strength B. In the thermo- 
dynamical limit N ^ 00 the model ([2]) has a quantum 
critical point at ft, = B/2J = he = I, and for < 7 < 1 
it belongs the quantum-Ising universality class [sit . 

Entanglement properties of the model have been stud- 
ied, particularly close to its phase transition 0. Atten- 
tion has been focused in the entanglement between two 
spins or the entanglement between a block of spins and 
the rest of the lattice. Spin-spin entanglement was stud- 
ied by means of concurrence and it was shown that a spin 
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is directly entangled to its neighbor spins and the range 
of spin-spin entanglement R (i.e. the maximal distance 
at which two spins are still entangled) depends on h and 
7 0) • In particular for the ground state R diverges 
at the factorizing field Q hf = \/\ — ^'^ [s^] while it is 
non universal and short ranged at critical point [3] • Block 
entropy was studied in terms of Von Neumann entropy 
and it was shown that it saturates as the size of the block 
increases for a non critical system, while it diverges loga- 
rithmically at the critical point with a universal behavior 
ruled by its conformal symmetry [3^. 

We consider the system defined in Eq.([2]) in the ther- 
modynamical limit and we focus our attention in a sub- 
system made of three spins in different positions in the 
lattice (not necessarily a block). We are interested in 
the entanglement shared between the three spins, and 
how such ME is distributed along the chain. Thus, as 
sketched in the previous section, we compute the Neg- 
ativity between a spin and the other two with respect 
to all possible bipartitions of the three spins. A similar 
analysis was done in [l8l | for three fermions out of a Fermi 
gas. 

THREE-SPIN ENTANGLEMENT IN THE XY 
MODEL 

Due to the complete integrability of the XY model, 
it is possible to calculate the reduced density matrix of 
any number of spins (see Appendix for the details). The 
density matrix is a function of the parameters, h and 
7, the temperature T and of the distances between the 
spins considered: a — j —i and f3 — k — j. We first study 
the case T = focusing on the range of the three-spin 
entanglement and then we consider the case of a block of 
three spins at finite temperature. 

Ground State Entanglement 

We first focus on the range of three-particle entangle- 
ment in the ground state and check whether spins distant 
enough such that they do not share spin/spin entangle- 
ment are nonetheless globally entangled. With this aim, 
we consider a spin and a block of two spins more distant 
than the spin/spin entanglement range (see Fig[T]a)) and 
we study the spin/block entanglement as measured by 
the Negativity. We find that the block of two spins may 
be entangled with the external spin, despite the latter 
is not entangled directly with any of the two spins sep- 
arately (see Fig. [2] upper panel). Hence the range of 
such spin/block entanglement may extend further than 
the spin/spin entanglement range. 

An interesting case is when there is no two-particle 
entanglement between of any the three spins, but still 
they share ME. Specifically, we consider the symmetric 




b) • • T t tT^ t • ■ 

FIG. 1: Configurations of spins described in the text and 
whose entanglement properties are presented in Fig. [5] We 
fix the range of spin/spin entanglement R < 3 such that spins 
at a distance d > 4 are not directly entangled, a) 'Clustering' 
two spins increases the range of entanglement (the scheme is 
symmetric also for spins on the left of the marked one), b) 
Symmetric configuration of spins such that no two-particle 
entanglement is present, but still the spins share multiparticle 
entanglement. 



configuration of spins drawn in FiglT]3. As shown in Fig[2] 
(lower panel), the Negativity between the central spin 
and the other two Afcentr and between an external spin 
and the other two Afsxt may be non-zero, despite the 
three spins do not share spin-spin entanglement. 

We shall see that spins in the configuration of Fig. 
[T] b) can share bound entanglement. To prove it, the 
idea is to resort the "incomplete separability" condition 
described in the first section. In fact from Fig. [5] we see 
that Mexi may be zero even if Mcentr is non-zero. Thus 
in such case the density matrix of the spins is PPT for 
the two symmetric bipartitions of one external spin vs 
the other two (t | TT ^^^d TT I T) and Negative Partial 
Transpose (NPT) for the partition of the central spin vs 
the other two. We remark that PPT does not ensure the 
separability of the two partitions, nevertheless the state 
is bound entangled. In fact if we could be able to distill 
a maximally entangled state between two spins then one 
of two previous PPT partitions would be NPT and this 
cannot occur since PPT is invariant under LOCC (20l.l33|. 

The scenario described above quantitatively varies if 
different values of 7 are considered. Since the range of 
spin-spin entanglement depends on 7 [s^] the distance 
between the spins at which ME exists and its range 
are also 7-dependent. For both the configurations we 
studied we found that ME is short ranged. However, we 
notice that such range diverges at the factorizing field h / 
(analogously to what occurs for spin/spin entanglement 
[37|). Remarkably, this holds also for BE. In fact we 
observe that the range of Mcentr is always greater than 
the range of Nexi- Hence approaching hf there are 
always configurations of spin far apart enough to share 
BE and for ft, ^ /i/ its range diverges. 

In summary if the spins are far enough to loose all 
spin-spin entanglement, some ME may be still present. 
The nature of the entanglement may change from free to 
bound entanglement by tuning the magnetic field h. 
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FIG. 2: T — Negativities between one spin and the other 
two Vs magnetic field h are shown for both configurations of 
Fig [1] We consider 7 = 0.5. In this case outside the inter- 
val marked by the two solid vertical lines ( 0.83 ^ h < 0.91) 
the range of spin-spin entanglement is R < 3 (for values of h 
inside this interval R grows due to the its divergence at fac- 
torizing field hf = \/l — 7^ ~ 0.86 [s^l)- For configuration a) 
(upper panel), Afsiock signals the spin/block entanglement for 
a distance d — 4. For values of h outside the vertical lines, the 
Negativity signals genuine spin/block entanglement. For con- 
figuration b) (lower panel) , both Msxt (solid line) and Mcentr 
(dashed line) are plotted. For values of h outside the vertical 
lines the three spins share no spin/spin entanglement, hence 
for non zero Mexi and Ncentr free multiparticle entanglement 
is present. The latter turns in to bound entanglement for val- 
ues of h such that only Mcentr 7^ and Msxt ~ (both on 
the left and on the right of the sohd lines). 



Thermal Bound Entanglement 

Quantum states must be 'mixed enough' to be bound 
entangled. In fact in a geometrical picture they may be 
located in a region separating free entangled states from 
the separable ones [40|. In our case, a source mixing is 
the trace over the other spins of the chain: as pointed 
out in the previous section, such loss of information may 
be enough to induce bound entanglement between the 
three spins. However if the spins are near enough the 
reduced entanglement is free. Thus in this situation it is 
intriguing to investigate whether the effect of the ther- 
mal mixing can drive the T = free entanglement to BE. 
To answer this question we consider a block of three ad- 



jacent spins and study the entanglement between them 
as a function of T. We see (Fig. |3]) that the entangle- 
ment between the two external spins is the most fragile 
increasing T, dying first at a certain Tc2', then at a higher 
temperature Tci the central spin looses its entanglement 
with each of the external ones. Thus for T > Tqi there is 
no two-particle entanglement between the spins, but still 
some ME is present, as detected by the Negativity Nexi 
and Mcentr- Eventually the Negativity vanishes at T^^^j 



and Ta/p 



In the region T^^^^ < T < Tj^^ 



the 



condition of 'incomplete separability' occurs: the density 
matrix is PPT with respect to the external bipartitions 
and NPT with respect to the central one and thus the 
spins share BE. 

We want to remark here that the failure of Negativ- 
ity to detect some entangled states does not invalidate 
the scheme described above. In fact, when all the Nega- 
tivities are eventually zero, the three-spin state may be 
still entangled, but such entanglement would be surely 
bound (since the state has positive partial transpose). 
This means that the region where bound entanglement 
is present might be wider than the one marked in Fig. 
[3l but, any case, such thermal bound entanglement will 
always separate the high temperature separable region 
from the low temperature free entangled one. 

This behavior shown for the Ising model is also found 
for the entire class of Hamiltonians In fact although 
the temperature at which the different types of entan- 
glement are suppressed is 7 dependent the qualitative 
behavior of the phenomenon is a general feature of the 
model 
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FIG. 3: Entanglement shared in a block of the three adjacent 
spins. We consider 7 = 1. In this case only nearest neighbor 
and next nearest neighbor spins are entangled (hence R — 2) 
at T = 0, 0] ■ The lines in the T — h plane indicate the tem- 
peratures at which the corresponding type of entanglement 
disappears. In the marked region ZVe^t < T < IVccret,- BE 
is present. 
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CONCLUSIONS 

We studied the entanglement shared between three 
spins in an infinite chain described by the anisotropic 
XY model in a transverse field. We analyzed the Neg- 
ativity between groups of spins as in Figll] demonstrat- 
ing that the corresponding block/block entanglement ex- 
tends along the chain over a longer range (Fig. ^ com- 
pared to the spin/spin case. Such type of entanglement 
persists at higher temperature (see Fig. [3]). 

In a recent paper, thermal BE was found for system 
made of few spins |41| . Here we proved the existence of 
bound entanglement shared between the particles, even 
in a macroscopic system. It occurs naturally in certain 
region of the phase diagram for the 'last' entangled states 
before the complete separability is reached. In this sense 
the BE bridges between quantum and classical correla- 
tions. Being the BE a form of demoted entanglement, 
we found that it appears when quantum correlations get 
weaker. At zero temperature bound entanglement ap- 
pears (see Fig. [2]) when the spins are sufficiently distant 
each others and as in the case of the spin/spin entangle- 
ment, it can be arbitrary long ranged near the factorizing 
field. At r 7^ BE is present when the system is driven 
toward a completely mixed state as the temperature is 
raised(Fig. [3]). It would be an interesting to further 
study bound entanglement in this context (f.i. focusing 
on other kinds of bound entanglement [i^). 

A more refined classification of ME could be done fol- 
lowing the scheme of Ref. fis'l , discerning different regions 
in the phase diagram in terms of different entanglement 
classes (GHZ or W). In Ref. 43| certain witness opera- 



tors able to distinguish some W- and GHZ- mixed states 
were discussed. We point out, however, that in our case 
such witness are not able to detect tripartite entangle- 
ment, not even if suitably generalized as Instead, 
the three tangle [j], that in principle could identify the 
GHZ class, is difficult to handle for generic mixed state 
and requires a hard numerical effort (some progress is 
achieved for low rank density matrices [4j|). 

It is plausible that increasing the size of the subsys- 
tem considered will increase the range of the ME. For 
instance, the range of spin/block entanglement (see con- 
figuration of Fig. [1] a) will increase if we consider a 
larger block. Hence, a single spin can be entangled with 
more distant partners, if one allow to cluster them into a 
large enough block. It would be intriguing to study how 
spin/block entanglement and, in general, block/block en- 
tanglement between subsystems, scale increasing the size 
of blocks, especially exploring the connection with quan- 
tum criticality. We remark that such analysis would be 
different with respect to the well known block entropy 
settingjl, li^, since in that case one is interested in the 
block/rest-of-the-system entanglement. Along this line, 
some results for a chain of coupled harmonic oscillators 



were obtained[46|. 

Finally, we speculate about the effect of temperature 
for the entanglement of a larger size block of spins. The 
behavior depicted in Fig. [3] suggests that entanglement 
between sub-blocks of few spins will be suppressed by in- 
creasing the temperature, before of the entanglement be- 
tween larger sub-blocks. Namely, temperature will sup- 
press first spin/spin entanglement, then spin/2-spinblock 
entanglement,.., n-spinblock/m-spinblock entanglement, 
and so on. Thus, at high enough temperature, parti- 
cles of the system will eventually loose all entanglement 
passing through successive steps at which temperature 
suppresses entanglement with respect to a 'microscopic' 
to 'macroscopic' hierarchy. 

We thank J. Siewert for significant help. We acknowl- 
edge fruitful discussions with A. Osterloh, F. Plastina. 
The work has been supported by PRIN-MIUR and it has 
been performed within the "Quantum Information" re- 
search program of Centro di Ricerca Matematica "Ennio 
De Giorgi" of Scuola Normale Superiore. 



Appendix 

It is convenient to express pijk in terms of three-point 
correlation functions that can be obtained explicitly fol- 
lowing the method used for the two-point correlation 
functions in 13611 



Pijk 



(3) 



In the previous equation p,q,r= 0, x, y, z (cr" = 1). 

Three-spin reduced density matrix is obtained from 
Eq.([3]). Due to the parity symmetry of the Hamiltonian 
47| . for the non-broken symmetry case some of the cor- 
relators are identically zero and the 8x8 matrix reads: 



Pijk 



\f Pi Phi 
8 V Pin Pii 



where 
Pi = 

/1 + A + B+ 





PII = 

(l+C-B- 



l + C + B- + 










1-C-B-+ j 



1-C-S+- 







V--£j^ l-C + B 






1-A + B++ J 



6 



Phi 



V- 



\ 



a/3 

V- 



Kp 





T^t Kp 







The entries of the matrices above are hnear combination 
of the correlation functions: 



A = 

±± _ 



B 

C = 

where f.i. YY ZaB — 



ZZ + ZZZ 

ZZd i ZZp i zz,^ 

Z -ZZZ 

XXZ^fj ± YYZ^p 
XZX^p ± YZY^p 



with a = j — i and (3 = 
k—j (the correlators depend only on the relative distance 
between the spins because of the translational symmetry 
of the system). 

Such three-point correlation functions can be calcu- 
lated following the method used in [s^ for the two-point 
ones. In short, after the Jordan- Wigner transformation 
which maps spin into spinless fermions 



CT" = - 



s=l 

l-l 

Bi W A,Bs 



s=l 

-AiBi . 



where A = 



") and B = 



the three- 



point correlation functions {crfcrja^) can be written as 
a Pfaffian whose elements are the two-point correlators 
(AiBm)- The structure of the Pfaffians further simpli- 
fies in to a determinant because only terms {AiB„i) — 
Gm-i = /J" ^{cos{(l)\m - ?|)(cos(0) - a) - ^s\n{(j)\m - 

tanh( ^ -y/ 7^ sin (l)'^~\-{h—cos 0-^)) 



i\)sH4>)}- 

ing. 



■\J 7^ sin (/)-^4-(/i — cos cfy^ ) 



are non vanish- 





Go Ga 


Ga+0 


ZZZafi — — 


G-Q, Go 


Gp 




G-a-f3 G-p 


Go 



XXZ^p = X 

G_i Go 

G_2 G_i 

G-Q G_Q_i 

G-a-p G-a-P-1 



Gq-2 Ga+P-1 
Ga-3 Ga+p~2 



G_i 

G-/3-1 



Gp 
Go 



XZXap = {-r+l'y. 

G—i . . . Ga-2 Go 



G, 



G-a+l 
G-a-1 



Go 
G_2 



G2 

Go 



a+p~2 

Gp 

G/3_2 



G-a-p ■ ■ ■ G-p-i G-p+i . . . G_i 



YYZ, 



aP 



Gi 
Go 



G2 
Gi 



G^a+2 G^a+3 
G-a-p+1 G 

~a — P+2 



Gq Ga+p 
Ga-1 Ga+P-1 

Gl Gp+i 

G^p Go 



YZY 



aP 



(-) 

Gl 



G_Q+2 

G-a 



a+P 



Ga-1 Ga 



+ 1 



Go 

G-2 



G2 
Go 



G-a- 



a-P+2 



G—p G- 



p l^-P+2 



Ga+P 

Gp+1 
Gp-i 

Gl 



The main difference with respect the two-point corre- 
lators calculated in [36] is that in this case the Toeplitz 
structure of the matrix is no longer valid. 
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